Tracing the evolution of the symmetry energy of hot nuclear fragments from the 
compound nucleus towards multifragmentation 
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The evolution of the symmetry energy coefficient of the binding energy of hot fragments with in- 
creasing excitation is explored in multifragmentation processes following heavy-ion collisions below 
the Fermi energy. In this work, high-resolution mass spectrometric data on isotopic distributions 
of projectile-like fragments are systematically compared to calculations involving the Statistical 
Multifragmentation Model (SMM). Within the SMM picture, the present study suggests a gradual 
decrease of the symmetry energy coefficient of the hot primary fragments from 25 MeV at the com- 
pound nucleus regime towards 15 MeV in the multifragmentation regime. The isotopic distributions 
of the hot primary fragments are found to be very wide and extend towards the neutron drip-line. 
These findings are expected to have important implications to the modeling of the composition and 
the evolution of hot and dense astrophysical environments, such as those of core-collapse supernova. 



PACS numbers: 25.70.-z, 25.70.Hi,25.70.Lm,26.30.+k 



Nuclear multifragmentation is one of the most interest- 
ing phenomena in nuclear physics as it holds promise for 
understanding nuclear matter properties at the extreme 
conditions of high excitation energy and large isospin 
(N/Z) asymmetry d 0, 0, EL H • The latter > in particular, 
plays a profound role in the dynamics of various astro- 
physical environments @, 0, 0, El ■ It is well established 
(e.g. [4[) that nuclear systems with relatively low excita- 
tion energy (E* j A < 2 MeV) form the traditional com- 
pound nucleus, whereas at higher excitation energy, the 
hot nuclear system expands and, subsequently, disassem- 
bles into an ensemble of hot primary fragments. This ex- 
tremely complicated process, namely the multifragmen- 
tation, occurs on a timcscalc of 100 fm/c (3.3xl0 -22 sec) 
during which the system can sample a large number of 
configurations. For this reason, statistical calculations 
(e.g., 0, [ll|]) have been very successful in describing 
this process. 

Recently, a remarkable similarity has been pointed 
out between the thermodynamic conditions (tempera- 
ture, density, isospin asymmetry N/Z) reached in nu- 
clear multifragmentation and the collapse/explosion of 
massive stars [ijL EH, Hi- This observation opens up 
the possibility of applying well-established models of nu- 
clear reactions to describe matter distribution and evo- 
lution during supernova explosions [HJ • In addition, sta- 
tistical calculations suggest that in multifragmentation 
Tl , an d i n hot astrophysical environments (e.g. su- 
pernova) [I2L [lH , the ensemble of primary fragments in- 
cludes neutron-rich nuclei towards or beyond the neutron 
drip- line. 

The primary fragments are expected to be hot (with 
excitation energies approaching 2-3 MeV/nucleon (l7j ) 
and, initially, in close proximity to neighboring frag- 
ments or nucleons. These conditions render their prop- 
erties, e.g. binding energy, different from those of cold 



(ground state) isolated nuclei. In particular, recent stud- 
ies [l8 , Til 20, 21 1 give evidence for a significant decrease 
of the symmetry energy of hot primary fragments. We 
recall that for a system with TV neutrons and Z pro- 
tons, the symmetry energy can be expressed as £ sym = 
C sym (N — Z) 2 /A, with C sym the symmetry energy coef- 
ficient and A the mass number. In the aforementioned 
studies, the symmetry energy of hot fragments from mul- 
tifragmentation at high excitation energy (E* /A = 4-6 

15 MeV or lower, 
25 MeV 



MeV) is found to be reduced to C sy 
as compared to the conventional value of C t 



sym 

for cold nuclei. In the same vein, our recent work [22| in- 
dicated decreased symmetry energy in the region E* /A 
= 2-3 MeV. The observed reduction of the symmetry en- 
ergy implies that more neutron-rich nuclei are favored in 
the distribution of fragments after the partitioning of the 
initial hot and expanded nuclear system. Such unusual 
neutron-rich nuclei should also be quite abundant in su- 
pernova matter [12j and may affect the dynamics of the 
explosion and the subsequent nucleosynthesis. 

For completeness, we emphasize that presently there is 
strong and growing interest in the density dependence of 
the symmetry energy of nuclear matter (e.g. P,0,[H, 24 1 
and references therein). However, we should point out 
that at low densities (less than half of the normal nu- 
clear density) the uniform nuclear matter breaks-up into 
fragments. Thus, in order to describe hot fragment dis- 
tributions, we should address the problem of the symme- 
try energy of individual fragments (primary fragments) 
embedded in a hot and dense environment. 

The present Rapid Communication is intended to pro- 
vide evidence for a continuous evolution of the symme- 
try energy coefficient from the compound nucleus regime 
(E*/A < 2 MeV) towards multifragmentation (E*/A> 4 
MeV) by a systematic comparison of the isotopic distri- 
butions of heavy fragments with calculations based on the 
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Statistical Multifragmentation Model (SMM). Herein, 
fragments with A > 40 are referred to as heavy residues, 
those with A < 40 as heavy IMF (intermediate mass frag- 
ments) , whereas both groups arc collectively called heavy 
fragments. Moreover, in this work, the N/Z distributions 
of hot nuclei from the multifragmentation of neutron-rich 
systems are found to be very wide involving exotic nuclei 
towards the neutron drip-line. 

The experimental data were obtained at the Cyclotron 
Institute of Texas A&M University employing beams 
from the K500 Superconducting Cyclotron. Two mag- 
netic separators were used: the Momentum Achromat 
Recoil Separator (MARS) for the 86 Kr (25 MeV/nucleon) 
+ 64 Ni, 124 Sn, 208 Pb reactions and the Superconducting 
Solenoid Line (BigSol) for the 64 Ni (25 MeV/nucleon) + 
64 Ni reaction. A detailed description of the measure- 
ments is presented in [i^, 25, [26}. Employing standard 



Bp-TOF-AE-E (magnetic rigidity, time-of- flight, en- 
ergy loss and residual energy, respectively) techniques 
[261 ] . the atomic number Z, the mass number A and the 
velocity of the fragments were obtained with high reso- 
lution. The measurements of 86 Kr + 64 Ni, 124 Sn were 
performed in the angular range 2.7°-5.4° and Bp range 
1.3-2.0 Tm. The measurements of the 86 Kr + 208 Pb were 
obtained in the angular range 1.0°— 2.7° and Bp range 
1.2-1.7 T m. Finally, the 64 Ni + 64 Ni measurements were 
performed with BigSol in the angular range 1.5°-3.0° and 
Bp range 1.1-1.6 Tm 22j. The measured yield distribu- 
tions were used to extract the average Z/A values. To 
obtain total fragment cross sections, the measured yield 
data were corrected for the spectrometer acceptance (as 
in [2^1) with the aid of the simulations described below. 

The calculations are based on a two-stage Monte Carlo 
approach. The dynamical stage of the collision was 
described by the deep-inelastic transfer (DIT) code of 
Tassan-Got 27J simulating stochastic nucleon exchange 
in peripheral and semiperipheral collisions. This model 
has been successful in describing the isospin, excitation 
energy and kinematical properties of excited quasiprojec- 

tiles in a variety of recent studies US H, Ha HI . The 

deexcitation of the quasiprojectiles was performed with 
the latest version of the Statistical Multifragmentation 
Model (SMM) 0, 0, EH , referred to as "SMM05" and 
briefly summarized below. This model assumes statisti- 
cal equilibrium at a low-density freeze-out stage and in- 
cludes all breakup channels ranging from the compound 
nucleus to vaporization (channels with only light parti- 
cles A<A) , allowing a unified description of nuclear disin- 
tegration with increasing excitation. In the microcanon- 
ical treatment, the statistical weight of a decay chan- 
nel is calculated as exponential of the entropy. Light 
fragments with A<A are considered as stable particles 
(nuclear gas) with only translational degrees of freedom. 
Fragments with A>4 are treated as heated liquid drops 
with free energies parameterized as a sum of volume, sur- 
face, Coulomb and symmetry energy terms 14| with pa- 



rameters adopted from the Bethe-Weizsacker mass for- 
mula. The model generates a Markov chain of parti- 
tions (with the Metropolis algorithm) representative of 
the whole partition ensemble [3l| . In the freeze-out con- 
figuration, the hot primary fragments are assumed to 
be isolated and at normal density; they, subsequently, 
propagate in their mutual Coulomb field, while undergo- 
ing binary deexcitation via evaporation, fission or Fermi 
breakup [lo| . The Coulomb interaction energy is di- 
rectly calculated for each spatial configuration of frag- 
ments in the freeze-out volume. Finally, the effect of the 
Coulomb field of the target in proximity to the decaying 
quasiprojectile is included. The symmetry energy evolu- 
tion is taken into account in the mass calculation of the 
primary fragments and the secondary fragments during 
the deexcitation stage, as in [L4| . Below a threshold of 
E* I A=\ MeV, a smooth transition to standard experi- 
mental masses is assumed [3]. Standard parameters of 
the SMM code are employed. A multiplicity dependent 
parametrization of the free volume [10j (determining the 
contribution of the fragment translational motion to the 
partition probability [l0() is used, whereas, the freeze-out 
volume (determining the Coulomb energy of the fragment 
partition) is taken to be 6 times the nuclear volume at 
normal density (as suggested by the asymptotic velocities 
of the observed heavy fragments). 

Fig. 1 shows the average proton fraction (Z/A) of the 
fragments with respect to mass for the reactions pre- 
sented in order of increasing system N/Z. The exper- 
imental data are shown by solid symbols and span the 
whole range of masses from evaporation residues down 
to heavy IMF. Vertical arrows indicate an approximate 
separation between heavy residues from evaporation and 
from the onset of multifragmentation. The thin solid line 
(marked "SL" ) gives the line of f3 stability and the thin 
dotted line (marked "EAL") represents the evaporation 
attractor line [12] (corresponding to the locus of fragment 
yields produced by evaporation of neutrons and charged 
particles from highly excited nuclei close to stability). 
The evaporation residue regime is populated mostly by 
neutron deficient products, except for the heaviest frag- 
ments from 86 Kr+ 64 Ni and 86 Kr+ 124 Sn corresponding 
to very peripheral products [25| . Fragments on the left 
side of the arrow (produced above the multifragmenta- 
tion threshold of E*/A ~ 2 MeV) show progressively 
smaller values of Z/A (Fig. 1) and cross the line of f3 
stability. 

The results of the DIT/SMM05 calculations are pre- 
sented by thick lines. The dotted lines correspond to 
calculations with constant values of the symmetry en- 
ergy of C sym = 25, 20, 15 MeV from the upper to the 
lower line, respectively. We observe that the standard 
calculation with C sym = 25 MeV produces fragments 
that are, on average, more neutron deficient than the ob- 
served fragments. With C sym = 20 MeV, the agreement 
with the data seems to improve, except for the lowest 
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FIG. 1: (Color Online) Average Z/A vs A. (a) 
64 Ni(25MeV/nucleon)+ 64 Ni (BigSol data). (b), (c), (d) 
86 Kr (25 MeV/nucleon) + 64 Ni, 124 Sn, 208 Pb, respectively 
(MARS data). Solid points: data. Thin solid line (SL): line of 
stability. Thin dotted line (EAL): evaporation attractor line. 
Thick lines: DIT/SMM05 calculations. Dotted lines, top to 
bottom: with C sy m = 25, 20, 15 MeV, respectively. Thick 
solid line: with C sy m{E* / A) given by the thick line in Fig. 2. 
Arrows: multifragmentation vs compound nucleus regime. 



masses, whereas with C sym = 15 MeV, only the latter 
group of masses is described. We point out that the 
measured heavy fragment data at 25 MeV/nucleon span 
a continuous range of excitation energies up to E* /A = 
4.0-4.5 MeV (corresponding to ~ 95% of all produced 
fragments). Assuming a continuous behavior of the sym- 
metry energy from low to high excitation, we tested var- 
ious forms of the dependence of C sym on excitation en- 
ergy E* I A. We found that the form represented by the 



E* I A > 4 MeV] provides a remarkable agreement of the 
calculations with the data from multifragmentation-likc 
processes, as shown by the solid lines in Fig. 1. 

The thin solid lines and the thin dashed line in Fig. 

2 indicate three additional forms of C sym (E* / A) that we 
tested. The results of (Z/A) for Kr+Sn are shown in Fig. 

3 by the corresponding thin lines (the lower one corre- 
sponding to the left thin line of Fig. 2 and the upper one 
to the right thin line). In Fig. 3, we notice that the thin 
solid lines embrace the thick line obtained with the "2-4" 
form of C sym (E* /A), whereas the thin dashed line is al- 
most indistinguishable from the thick one. Similar results 
are obtained for the other reactions of this work. Fur- 
thermore, assuming that the multifragmentation thresh- 
old is at E* j A = 2 MeV, we may only adopt forms of 
C sym (E* I A) with the condition C sym = 25 MeV at the 
threshold and below. An appropriate such form appears 
to be the "2-4" form represented by the thick solid line 
in Fig. 2, as mentioned previously. Consequently, within 
the present SMM interpretation of fragment formation, 
the symmetry energy coefficient appears to decrease from 
25 MeV in the compound nucleus regime (left arrow in 
Fig. 2) to 15 MeV in the bulk multifragmentation regime 
(right arrow in Fig. 2) . This phenomenological evolution 
of the symmetry energy of hot heavy fragments is the 
main result of the present study. 

For comparison, in Fig. 2 we include the values of 
the symmetry energy coefficient obtained in our recent 
work on heavy fragment isoscaling from the 86 Kr (cir- 
cle, square) and 64 Ni (triangle) reactions We notice 
that, in spite of the different approach of symmetry en- 
ergy determination, an overall agreement in the trend is 
observed. We speculate that the reduction of the sym- 
metry energy with excitation observed in this work and 

the aforementioned previous studies [3 m, m mi may 

originate from in-medium modifications of the properties 
of the hot primary fragments in their dense surrounding 
Detailed microscopic calculations are neces sary 



to shed light on this interesting and challenging issue [34 1 . 

Apart from the average Z/A properties, the present 
calculations describe successfully the characteristics of 
the observed fragments including the widths of the iso- 
topic mass distributions, the velocity distributions and 
the production cross sections. A full report on these re- 
sults is currently underway. For completeness, we present 
in Fig. 4 a comparison of the experimental mass distri- 
butions for several elements with the DIT/SMM05 cal- 
culations for two reactions: 86 Kr+ 64 Ni and 86 Kr+ 208 Pb. 
The "2-4" form of C sym (E*/A) dependence (thick line 
in Fig. 2) was used in these calculations. For Kr+Ni, 
the data are shown by the full circles and the calcula- 
tions by the thick solid lines (the thin solid lines giv- 
ing the hot primary fragment distributions for these el- 
ements). Similarly, for Kr+Pb, the data are given by 
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FIG. 2: (Color online) Thick solid line: form "2-4" of 
C syTa (E* /A) in the SMM05 calculations. Thin lines: other 
forms of C ay m{E*/A) tested (see text and Fig. 3). Thin 
horizontal lines (top to bottom): C sym = 25, 20, 15 MeV, 
respectively, used in Fig. 1 calculations. Solid points: values 
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FIG. 3: (Color Online) (Z/A) vs A for 86 Kr + 124 Sn. Solid 
points (data), lines SL, EAL and arrow as in Fig. 1. Thin lines 
tracing the data: DIT/SMM05 calculations with C sym _(E*/A) 
forms presented by the respective lines in Fig. 2 (see text). 
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FIG. 4: (Color online) Comparison of experimental mass 
distributions with DIT/SMM05 calculations. Solid points: 
86 Kr+ 64 Ni data. Open points: 86 Kr+ 208 Pb data. Full lines: 
calculations for 86 Kr+ 64 Ni; thick: final (cold) fragments, 
thin: primary (hot) fragments. Dotted lines: calculations 
for 86 Kr+ 208 Pb; thick: final (cold) fragments, thin: primary 
(hot) fragments. Arrows: neutron-drip line [35| . 



open circles and the calculations by thick dotted lines 
(again, the thin dotted lines describing hot fragments). 
We find an overall satisfactory agreement of the calcula- 
tions with the mass distribution data. Interestingly, the 
Kr+Pb data extend to more neutron-rich products than 
those of Kr+Ni (demonstrating the sig nificant role of the 
target N/Z at this energy regime [25|, |3fJ), a trend well 
described by the calculation. Finally, we point out that 
the distributions of the hot primary fragments are wide 
and include very neutron-rich nuclei. Particularly, for 
the Kr+Pb reaction (the most neutron-rich system stud- 
ied) , they are very broad and extend towards the neutron 
drip-line (35| (arrows in Fig. 4). Our calculations indi- 
cate that the distributions of these hot exotic nuclei are 
sensitive to their masses and, in turn, sensitive to the 
evolution of the symmetry energy with excitation. 

We wish to emphasize that high-resolution mass spec- 
trometric data on heavy fragment (heavy residue and 
heavy IMF) distributions - spanning the full spectrum of 



reaction mechanisms (and excitation energies) from the 
traditional compound nucleus towards multifragmenta- 
tion - may provide a promising experimental probe of the 
properties of hot and extremely neutron-rich fragments 
when coupled with detailed calculations as presented in 
this study. From an experimental point of view, we wish 
to propose systematic measurements on residue distribu- 
tions on a variety of neutron-rich systems. The energy 
regime of ~25 MeV/nucleon has two main advantages: 
(a) N/Z equilibration is attained (e.g. 30]), thus the 
target has maximum contribution to the N/Z distribu- 
tion of the quasiprojectiles and (b) the excitation energy 
is not too high (E* /A < 4 MeV), so that hot heavy pri- 
mary fragments are abundantly present. A very exciting 
extension of these studies in the future involves reactions 
with very neutron-rich rare beams on heavy neutron-rich 
targets in current or planned radioactive beam facilities 
(e.g. the Advanced Exotic Beam facility in the US). 
We believe that comparisons of high-resolution heavy 
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fragment data with precise model calculations can pro- 
vide valuable information on the distributions and the 
properties of hot and very neutron-rich nuclei towards 
and beyond the neutron drip-linc. Knowledge on the 
properties of exotic nuclei produced in hot and dense 
environments are essential for astrophysical calculations, 
such as the composition and dynamics of core-collapse 
supernova and the course of the relevant nucleosynthesis 
processes (e.g., the r process). 

In summary, experimental data on isotopic distribu- 
tions of projectile-like fragments from peripheral and 
semiperipheral collisions of 25 MeV/nucleon 86 Kr and 
64 Ni beams on heavy targets are compared to calcula- 
tions involving the Statistical Multifragmentation Model 
(SMM). Within the adopted SMM picture, the present 
study suggests a gradual decrease of the symmetry en- 
ergy coefficient of the hot primary fragments from 25 
MeV in the compound nucleus regime towards 15 MeV 
in the multifragmentation regime. The isotopic distribu- 
tions of the hot fragments following the multifragmenta- 
tion stage are very broad and, for the most neutron-rich 
systems, are found to extend to exotic neutron-rich nu- 
clei. These findings are of importance in calculations of 
the isotopic composition in hot and dense astrophysical 
environments, such as those encountered in supernova. 
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